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We study braneworlds in a five dimensional bulk, where cosmological expansion is mimicked by 
motion through AdSs . We show that the five dimensional graviton reduces to the four dimensional 
one in the late time approximation of such braneworlds. Inserting a fixed regulator brane far from 
the physical brane, we investigate quantum graviton production due to the motion of the brane. We 
show that the massive Kaluza-Klein modes decouple completely from the massless mode and they 
are not generated at all in the limit where the regulator brane position goes to infinity. In the low 
energy limit, the massless four dimensional graviton obeys the usual 4d equation and is therefore 
also not generated in a radiation-dominated universe. 
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In recent years, cosmological models where our Universe is represented as a hypersurface moving through a higher 
' dimensional spacetime have received a lot of attention. The reasons for that are twofold. First, string theory, the 
presently most promising candidate of a theory for quantum gravity is consistent only in more than 3 + 1 , namely in 
9 + 1 or (for M-theory) 10+1 dimensions G|, |J. Second, the hierarchy problem i.e. the unnatural disparity of the 
electroweak scale ~ ITeV and the Planck scale ~ 10 16 TcV can be addressed in the context of extra-dimensions ||. 
■ For this, it is necessary that one or several of the extra-dimensions be much larger than the Planck scale. The 4- 

dimensional effective Planck mass then becomes M 4 ~ M%+™L n , where n denotes the number of extra dimensions 
and L is their size. These extra dimensions can be large and nevertheless undetected if standard model particles are 
ly-j ' confined to a lower dimensional hypersurface and cannot probe them. String theory predicts the existence of such 
Dp-branes onto which standard model particles are confined Q. The bulk spacetime around the Dp-brane can then 
""^5 ■ only be probed by gravity. Gravity has been tested only down to scales of about 0.1mm and therefore allows for 
L < 0.1mm. It has been shown that confinement of gravity to a region of size L < 0.1mm around a brane with one 
extra-dimension can also be achieved by a non-compact Anti-de Sitter spacetime ||. The cosmological situation of 
I an expanding universe is then obtained by a brane moving through a 5-dimensional Anti-de Sitter spacetime ( AdS 5 ) . 
At low energy, this setup leads to the usual Friedmann equations for the expansion of the Universe || [?J . This is the 
case which we consider in this paper. 

We want to address the following issue: a brane moving through AdS 5 spacetime leads to time-varying boundary 
conditions for quantum fields living in the bulk. In quantum physics it is well known that moving boundaries yield 
particle creation from vacuum, so-called motion induced radiation For instance, photons are produced from 
vacuum in dynamical cavities (dynamical Casimir effect) H as well as by a single moving mirror ]l0| . Here we study 
the same mechanism in brane world cosmology, namely graviton generation provoked by the motion of the brane. 
Apart from the massless graviton, braneworlds allow for a tower of massive Kaluza-Klein gravitons which might also 
be produced by the moving brane. Such massive modes have potentially devastating effects as they would eventually 
dominate the energy density of the universe and spoil the phenomenology if their production is sufficiently copious. 
Here we shall show explicitly that the dangerous massive modes are not produced in the single brane case. Also the 
production of massless 4-dimensional gravitons is strongly suppressed at low energy. 

In the next section we present our setup and derive the equations of motion for the graviton modes. We espe- 
cially discuss the nature of the coupling matrix which describes the coupling of the different modes due to the time 
dependence of the boundary (the brane). The formalism laid out is very general and we plan to apply it to other 
setups (high energy, 2-brane models) in the future. We also show that at low energy and for a regulator brane placed 
sufficiently far, the massless graviton mode obeys precisely the 4d equation of motion for gravitons in a Friedmann 
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universe and decouples from the massive modes. In Section II] we brielfiy estimate the produced gravitons and their 
spectrum analytically, focusing on the dominant 0-mode. We summarize our conclusions in Section |v|. 

II. THE MODE EQUATIONS FOR GRAVITONS IN MOVING BRANEWORLDS 

A. The background 

We consider the cosmological Randall-Sundrum-II model where we have AdS 5 geometry in the bulk with one physical 
brane at a time-dependent position yb(t)- In Poincare coordinates, the bulk metric is given by 

J 2 

ds 2 = g AB dx A dx B = [-dt 2 + S ij dx i dx j + dy 2 ] . (2.1) 

Capital Latin indices A,B run from to 4 and lower case Latin indices i,j from 1 to 3. Four-dimensional indices 
running from to 3 will be denoted by lower case Greek letters. L is the AdS 5 curvature radius and is related to the 
negative cosmological constant A by A = — Q/L 2 to solve the bulk Einstein equations 

G AB + Ag AB = . (2.2) 

We now introduce a brane at y = yb(t) and replace the "left hand side", < y < yb{t), of AdS 5 by a second copy of 
the "right hand side". We use the superscripts ">" and "<" for the bulk sides with y > yb and y < —yb, respectively. 
In terms of the coordinate y, the value of y decreases continuously from oo to yb and then jumps to — yb over the 
brane whereafter it continues to decrease. At the brane position, y^ = yb(t), y^ = — j/b(^), the metric function (L/y) 2 
has a kink. The Einstein equations at the brane position are singular, they contain a Dirac-delta function, 

Gab + AgAB = K 5 T%T e 5(y - y h ) . (2.3) 

The delta function confines the e nerg y momentum tensor from Standard Model fields to the brane. To avoid the delta 



function, one can integrate Eq. (2.3) over the extra dimension which leads to the so-called Israel-Darmois junction 
conditions (ll], 12, [l3| [l4j at the brane position. These read |H| 

9% - 9tu = , (2.4) 

where S^ v is the energy-momentum tensor on the brane with trace S, and 

Kr ee 6n 2 G 5 = -L . (2.6) 
M 5 and G5 are the five-dimensional (fundamental) reduced Planck mass and Newton constant, respectively. K uv i s 



the extrinsic curvature of the brane and is the induced metric on the brane. The first junction condition (2.4) 
simply states that the induced metric, the first fundamental form, 

q^v = e^g AB , (2.7) 

be continuous across the brane. Here the vectors are tangent to the brane. If we parametrize the brane by 
coordinates (z M ) and its position in the bulk is given by functions AT^(z M ), the vectors e M can be defined by 

e£ = d,X£(z) . (2.8) 

Denoting the brane normal by n, we have g AB e^n B — 0. The extrinsic curvature can be expressed purely in terms of 
the internal brane coordinates [H], 17 1, K — K ^dz" 'dz^ 1 , with 



K„ v = -\ [g AB (e*d v n B + e^n B ) + e*eln c g AB , c \ . (2.9) 

The link between the extrinsic cu rvat ure of the hypersurface and the brane energy-momentum tensor is established 
by the second junction condition (|2.5|), which replaces the 4d Einstein equation. 
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A homogeneous and isotropic brane moving through AdS 5 with brane position j/b(f?) and bulk time given by t\,(j]) 
has the metric 



ds 2 



1 



^] \dl 2 +,\ i d.r'a.r> 



a 2 {ij) [-dry 2 + dx l dx J ] , 



where a = L/y^ is the scale factor and r\ denotes the "conformal time" of an observer on the brane, 



(2.10) 



(2.11) 



The brane motion induces the 7-factor which relates the (conformal) eigentime r\ of the brane to the (conformal) 
coordinate time t of the bulk. From now on an overdot indicates a derivative w.r.t. conformal time rj on the brane. 
The brane normal is then given by 



n 



-Hi 



tv. (ty. 



Vb 



(2.12) 



We consider a homogeneous and isotropic total energy momentum tensor on the brane, Stf = — TStf. Here T 
is the brane tension and T{? is the energy momentum tensor of particles and fields confined on the brane given by 



-p, Tj = PSj. The second junction conditions now become 



K 5 (p + T)=6 



K 



x/i + L 2 H 2 

L ' 
2LH 
ay/1 + L 2 H 2 
3Ha(p + P) , 
2 



H< = f 8 Tp[l 



2TJ 
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(2.13) 

(2.14) 
(2.15) 
(2.16) 



where H = a /a 2 . Equations ( |2.13|) to 
the literature (for a review, see |18| or 



( 2.16| ) form the basis of brane cosmology and have been discussed at length in 
19f|). The last equation is called the "modified Friedmann equation" for brane 
cosmology J/J. For usual matter with p + P > 0, p decreases during expansion and at sufficiently late time p . 
The ordinary 4-dimensional Friedmann equation is then recovered if we set 



k 2 T 2 



Neglecting the 4-dimensional cosmological constant, A4 ~ 0, we obtain in addition 

6 K r 

— — and k. = — . 

k,T 4 L 



and k 4 = 8ttG4 = 



n 2 T 



L = 



(2.17) 



(2.18) 



Note that, although for a de Sitter brane the density is simply constant, p = —P = est and there is no late time 
approximation, Eq. ( 2.14 ) implies that the Hubble rate remains constant and we reproduce the usual exponential 
expansion. Only the relation between the expansion rate H and the brane density p is modified. 



B. Tensor perturbations 



A quantum field in the bulk is generically expected to be modified by the moving brane which acts as a moving 
boundary of the 5-dimensional bulk spacetime. We want to study this effect for bulk gravitons. With this in mind, 
we now linearly perturb the bulk metric allowing for tensor modes 1 , 



ds 2 



I? 
1,2 



-dt 2 + {S i:j + 2h lJ )dx i dx j + dy 2 ] 



(2.19) 



1 The word "tensor" here has to be understood as spin 2 degree of freedom w.r.t. rotations of the 3-dimensional homogeneous and isotropic 
slices normal to the extra-dimension y. 
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Tensor modes satisfy the traceless and transverse conditions, h\ = = 0. We then decompose hij into spatial 
Fourier modes, 



hij(t,x,y) 



(2tt)3/2 ^ 



J2 e ih -e^h.(t,y,k) 



(2.20) 



where e'j are unitary constant transverse-traceless polarization tensors which form a base of the two polarization 
states . = + and . = x. The perturbed Einstein equations now yield the equation of motion for the mode function 
h., which obeys the Klein-Gordon equation for minimally coupled massless scalar fields in AdS 5 



21 



df + k 2 -d 2 + -d y 



h.{t,y;k) = 



(2.21) 



In addition to the bulk equation of motion the modes also satisfy a boundary condition at the brane coming from the 
second junction conditions, 



LHd t h. - a/1 + L 2 H 2 d y h. 



7 1 (vd t - d y ) h. 



(2.22) 



Here n^ T) denotes possible anisotropic stress perturbations in the brane energy-momentum tensor and we have intro- 
duced the brane velocity v = -7=====. For simplicity we set II^ T) = in this work. (Some of the difficulties which 

appear when I1^ T) ^ are discussed in [^3) . The wave equation ( [2.21 ) with boundary condition ( 2.22] ) cannot be solved 



analytically except if the background metric functions are separable, and this only happens for maximally symmetric 
branes, i.e., branes with constant Hubble rate. This includes the Randall- Sundrum case H = 0. A cosmologically 
relevant case is the de-Sitter brane, < H = est. The spectrum of gravitational waves generated during de Sitter 
brane inflation can be calculated [E2l E4L E5l E6| . 



C. Late time approximation 

We want to investigate late time and see whether we have graviton production also at late time, after inflation. We 
therefore consider the limit 

v = - LH « 1 . (2.23) 
Vl + L 2 H 2 

In this limit the boundary condition ( 2.22| ) reduces to 

d v h.\ yh =0. (2.24) 

If the position of the brane is fixed, the solutions of the system formed by Eq. ( 2.21| ) and Eq. ( [2.24 ) are well known. 
These are the Bessel functions, 

h. = Aexp (±iut) (my) 2 [J 2 (my) + BY 2 (my)} , (2.25) 



where ui = V m 2 + k 2 . The junction condition ( 2.24 ) requires 

R Ji(myb) 7T 2 

B = -Y[(^yV)^4 imyh) ' (2 - 26) 

where the last expression is a good approximation for my b -C 1. This is precisely the result of Randall and Sundrum j| 
for a static AdS 5 brane. Since the extra dimension y is not compact, the mass-spectrum is continuous, m 2 can take 
any non-negative value. 

Allowing for a single moving brane is not well suited for a numerical treatment. We therefore introduce a second, 
so-called "regulator brane" far away from the physical brane, at the position y — y T which we let tend to infinity 
at the end of our calculation. We assume the regular brane to be empty and fixed. The boundary condition at the 
regulator brane is thus 

d y h.\ y =0. (2.27) 
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For a solution of the form ( [2.25 ) this implies the additional constraint 

Ji(im/ r ) +BY 1 (my z ) = . 



(2.28) 



This condition is satisfied only for a discrete series of mass eigenvalues m a . With B also m a depends on the position 
of the physical brane and therefore on time. 

The evolution equation (2.21) together with the boundary conditions form a Sturm-Liouville problem, with eigen- 
value equation 



~dy + ~dy 

y y 



-y 3 d y [y 3 d y (j) a ] = m 2 a cj) a . 



At any given time, Eq. (2.29) yields an orthonormal system of "instantaneous solutions" 

(f>o=A Q + B Q y 4 , 

<j>i = A i (m l y) 2 [J 2 (m i y) + B^im^)} = Ai(m l y) 2 C 2 (Toa i y) , 



(2.29) 



(2.30) 
(2.31) 



where we have introduced C v {mj,y) = J u (miy) + BiY v {m.iy). Here and below we use indices i and j to denote 
all the mas sive modes and indices a and [3 to denote all the modes, including the massless mode. The boundary 
conditions (2.24) and ([2.27]) then require 



B Q = and C*iKt, b ) = d{m iyi .) = . (2.32) 
Furthermore, the solutions form an orthonormal system of functions on the Hilbert space J$? — Jz? 2 ([yb, Vr], 2/ _3 dy) 



(0a, 



V 3 (t> a (t,y)(l)[3{t,y)dy = 6 a p 



or, explicitly, 



AiAj 



< - {m J yC2(my)Ci(m j y) - m i j/Ci(m i y)C 2 (m :) ?/)]P' = 

m i — m. »° 

i,<h) = -AiiOi [jT 1 {Aadiuiiy) - B y A C 3 (miy)}] |^ = , 

'^^) = ^(^ m 2 [( m iy) 2 {Ci(miy) - Ci(miy)C 3 (miy)}] |^ = 1 , 



A 



2 o y ~h4 + AoBo(y 2 y 2 h ) + \b*{ v * -yl) = l 
z y r Vb D 



(2.33) 

(2.34) 
(2.35) 
(2.36) 

(2.37) 



The orthogonality relations are trivially satisfied with the boundary conditions ( 2.32| ), whereas the normalization 
conditions fix the constants A a , 



- (Atrai) 2 (m^y^Cf (imy) 



lift 



2 2 
A 2 Vi ~ Vb _ -, 
A ° 9„2„2 ~ 



^yivi 



(2.38) 
(2.39) 



All the constants depend over the brane position yh(t) on time. Since the orthonormal set {<f> a (t, y)}^L * s complete 
in Jif, we can expand a generic solution in the form 



h.{t,y;k) = ^q a {t]k)(j) a {t,y) 



q a (t;k) = (4> a ,h.) 



(2.40) 



a=0 



This enables us to write the second order action in h. leading to the equation of motion for h, as an action for the 
canonically normalized coefficients q Q . Denoting the second order perturbation of the gravitational Lagrangian by 
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5(y/gR) we obtain for the action 

= &J dtdyy ~ 3 ^ dth ' )2 ~ {dyH ' )2 ~ ^ 

= ^ J dt ^2 [( d t<la){dtqf3)Sai3 + 2M a pq a (d t qp) + {N a p - (dy(j) a ,d,, 

= \ f dt Yl [( 9 *q«) 2 - ^Wa + {2M a pq a d t qp + N a0 q a q } 

where we have introduced coupling matrices 

M a(} (t) = (d t <j) a ,<pp) i N a p(t) = (d t <p a ,d t (f>t 
and the canonically normalized coefficients 



k 2 S a p} q a qp] 



q a 



■ q a 



(2.41) 



(2.42) 



(2.43) 



As before, the frequency is given by 



tul{t-k) = k 2 + m 2 a 



(2.44) 



The equations of motion for the Fourier coefficients q a (t; k) then simply follow from the Euler-Lagrange equations for 
the action ( p. 41 ), 



d 2 t q a + L0 2 a q a + J2 i M p a - M a p] d t qp + [9 t M a0 - N a0 ] qp = 
P 



(2.45) 



Since ((j) a ,(j)p) = constant, one might naively expect that M a p be antisymmetric. But this is not the case, since 
M a p depends on ti me n ot only via <f> a and <fip, but also via the integration boundary which enters the definition of 
the inner product (2.33). On the other hand, the completeness of the eigenf unctions <j) a implies 



X! 4>i(y)4>i(y) = % - y)y 3 



(2.46) 



so that 

J^M^Mp-t = ^rXft^s.^y)] = {d t (t>a,d t (t> ) = N a p . (2.47) 

7 7 

The system of coupled second-order differential equations ( 2.45| ) thus solely depends on the time-dependent fre- 
quency u)\ and the coupling matrix M a p. Numerical solutions of this system in other situations (oscillating cavities) 
have been obtained in Refs. (2?|, [28[| . Also the results obtained in this work are tested with the code described in |p8| . 
We plan to apply this numerical technique in future work. 

The matrix elements are given in terms of the instantaneous solutions <p a . The use of several identities of Bessel 
functions finally leads to 



Moo = yb 9 Vl ' 2 

v 2 - vi 

M oj = , 



M i0 = 2Vm~Moo , 

Ma = ihj , 

Mij = .l,.l,m,m/mj / y 2 C 1 (m i y)C (m j y)dy 



2m? 



2 2 

mf — m ■ 

1 j 















l+2i" 




V 









! mj(y h + m 3 ) 
%(j/b + m) 



(2.48) 

(2.49) 

(2.50) 
(2.51) 



(2.52) 
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where we have set 



Vb(t) 



d t \n (mi) 
d t In (y h ) 



Vh 



Y 2 {u\,y h ) 



-Ha 



and 



-H 



(2.53) 
(2.54) 



Numerics indicate that both terms appearing in My have similar amplitudes, thus the second term can give us an 
order of magnitude estimate of Eq. (2.52). We are mainly interested in the limit y r — » oo. Setting e = yb/yr we find 
- 2 fore 



0. To lowest order in e and in the limit v <C 1 the matrix elements are given by 

Moo = -H[l + 0(e)] , 
M 0j = , 
M i0 = HO(e) , 
Ma = HO(e 2 ) , 
My = HO{e 2 ) . 



(2.55) 
(2.56) 
(2.57) 
(2.58) 
(2.59) 



This leads to the following interesting result. At late times t ~ rj and in the 1-brane limit e — > 0, the Kaluza-Klein 
modes with non-vanishing mass evolve trivially and only the massless mode is coupled to the brane motion 2 



qo 



k 2 -H-H 2 



qo = , 



q, + [k 2 + m 2 ] q 4 = , for i ^ . 
Therefore bulk tensor perturbations h. ~ h. Q = \f2{L j a)qo satisfy the equation 

h, + 2Hh. + k 2 h. = . 



(2.60) 
(2.61) 

(2.62) 



This equation is valid everywhere in the bulk, and in particular on the physical brane where it reproduces exactly 
the equation for a 4-dimensional gravity wave in a Friedmann universe. This explicitly proofs that at low energy 
(late times) the homogeneous tensor perturbation equation in brane cosmology reduces to the 4-dimensional tensor 
perturbation equation. The massive Kaluza-Klein modes decouple completely from the massless mode in the 1-brane 
limit. The friction term arises from the time-dependent boundary condition, hence from the motion of the brane, and 
is not due to the excitation of the massive modes as claimed in |29j. 

This is our first main result. In the cosmological 1-brane model (RSII with a moving brane), the 4-dimensional 
graviton obeys the usual equation of motion and the (canonical) massive modes behave like massive particles in 
Minkowski space. Contrary to the massless mode, they are not affected by the motion of the brane. However note 
that k and m are the co-moving momentum and mass. 



III. QUANTUM GRAVITON PRODUCTION 



So far, we have considered the classical problem of solving a partial differential equation with time-dependent 
boundary conditions. Now we go on to consider the quantum production of gravitons. Let us assume (for the sake 
of simplicity) that the physical brane yb has been at rest up to some initial time t- m ~ r) m ^ L/a- m , and that at this 
initial time there have been no gravitons. The graviton field has been in the vacuum state. Now the brane starts 
moving and (as seen from an observer on the brane) the universe expands. Most probably, at some later time ifi n , the 
graviton field will no longer be in the vacuum state, gravitons have been generated. 

Of course since the expansion of a radiation-dominated universe is decelerating already prior to ij n the universe 
is more likely to expand faster than to be at rest, hence our initial conditions arc not realistic. To study the 



quantum generation of gravitons we now consider the field h,(k). In addition to Eq. (2.4C), we have a similar mode 
decomposition for its canonical momentum, 



■K.(t,y;k) = \\— . = J — 5^Pa(*; k) <j> a (t,y) = p Q (f; k)<f> a (t, y) 

5 \ ') V 5 a=0 a=0 



(3.1) 



It is easy to see that e 1 ^2 a M a j is bounded for all values of e and therefore we do not have to fear that the infinite sum might 
contribute in the limit e — > 0. 
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The last equality defines the expansion coefficients p Q of the momentum n.. In quantizing the field, we promote both 
the field h. and its conjugate momentum tt. to operators. We then impose the equal time canonical commutation 
relations for bosonic fields, 



h.(t,y;k),h.(t,y';k) 
[tt. (t,y;k),%.(t,y';k)) 
h.(t,y;k),TT.(t,y';k) 



= , 
= , 

= iS(y-y') . 



The operator valued coefficients q a , p a thus satisfy 



[q a (t;k),qp(t;k)] = , 

[p<*(*;fc),p/3 = ' 

[qa(t;k),f>p(t-,k)] = iS a 



(3.2) 
(3.3) 
(3.4) 



(3-5) 
(3.6) 
(3.7) 



As usual in quantum field theory, we wo rk in the Heisenberg picture. The operator valued coefficients q Q (i; k) still 



satisfy the classical equation of motion ( 2.45 ) 
operators are given by 

q a (t;k) 



But now we also know the initial state. Before i; n the coefficient 



for t < ti 



(3.8) 



Here d Q (fc) and a^ik) are the annihilation and creation operators associated with the initial vacuum state. At time 
t, we define the vacuum state |0, t) as the state annihilated by all annihilation operators, a a (k, t) |0, t) = 0, V/c and a. 

We now consider a fixed wavenumber k, as before we work in the late time limit, hence r/ ~ t and we assume a 
barotropic fluid P = wp with w — est on the brane, so that the scale factor is a ~ (r)/L) 2 ^ 1+3 ' w \ We are mainly 
interested in a radiation-dominated universe, w — 1/3. In the limit where e — > , we may neglect the intcrmodc 
coupling. The massive Kaluza-Klein modes then obey the massive wave equation ( |2.61 ) with constant mass, hence 
they do not experience any particle generation. The massless mode evolves according to equation 



q + [k 2 - (v 2 - §) V - 2 ] q = 



3(1 - w) 
2(1 + 3u>) 



(3.9) 

\)n-* = H+H* = 



This is a standard harmonic oscillator equation with a time-dependent mass term m 2 = — (i 

a/a. In a radiation-dominated universe v = 1/2, the mass vanishes identically (the expansion velocity d is constant) 
and there is no particle creation. During a matter dominated universe m 2 (r]) is negative and time-dependent, so 
there is particle creation, but one easily estimates that the produced energy density is of the order Ph(j]) / Pc(ij) — 



l\l '(a V , 



cqj 



(h/r eq y 



10" 



which is vanishingly small. Here £^ — J~R~ A denotes the 4-dimensional Planck scale 



and T eq is the age of the universe at equality. 



IV. CONCLUSION 



We have shown that in the single brane cosmological Randall-Sundrum model at late time, no massive gravitons are 
produced due to the motion of the brane through AdS 5 . Even though we have only shown this for the 4-dimensional 
spin-2 mode, the same result is expected for the spin-1 "gravi-photon" and the spin-0 "gravi-scalar" modes which 
obey the same equations (they are the three additional helicity states of a massive spin-2 particle). Also the massless 
gravi-photon and gravi-scalar are not produced in the single brane limit since they become non-normalizable. The 
only particle which can be produced in this model is the massless spin-2 graviton. 

At late time, the massless spin-2 graviton obeys exactly the 4-dimensional equation of motion. Therefore, as in 
4-dimensional cosmology, there is no graviton production during the radiation-dominated phase of the universe. This 
is a very general result in 4-dimensional cosmology: during the radiation era, minimally coupled massless particles 
are not generated, however, if the universe obeys a different expansion law they are. 

Models with two branes which approach each other and later move apart again (e.g., the ekpyrotic model), might be 
more interesting for the study of graviton production. Due to mode coupling, we expect particle production also for a 
radiation-dominated expansion law in this case. Furthermore, in two brane models also massive modes are produced, 
which can lead to stringent constraints, as they soon come to dominate the universe. Unfortunately, the high energy 
regime, where most of the particle production takes place, depends on the details of the model under consideration. 
We will discuss the two brane case in a forthcoming paper |3C[|. 
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